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Abstract. In this paper we show how to obtain decay estimates for the 
damped wave equation on a compact manifold without geometric control via 
knowledge of the dynamics near the un-dampcd set. We show that if replacing 
the damping term with a higher-order complex absorbing potential gives an 
operator enjoying polynomial resolvent bounds on the real axis, then the "re- 
solvent" associated to our damped problem enjoys bounds of the same order. 
It is known that the necessary estimates with complex absorbing potential 
can also be obtained via gluing from estimates for corresponding non-compact 
models. 



1. Introduction 

On a compact, connected Riemannian manifold without boundary (X, g), we 
consider the non-selfadjoint Schrodinger operator 

(1.1) P(h) = h 2 A g +iha 

where a € C°°(X) is a non- negative function, and A g = d*d is the non-negative 
Laplacian associated to the metric g. This paper mainly addresses the question of 
the semiclassical analysis of the resolvent of P(h), 

R z (h) := (P^-z)- 1 

for z in a complex h— dependent neighborhood of 1. For non-selfadjoint operators, 
it is well known that the norm of the resolvent ||i? z (/i)||z;(i 2 ,L 2 ) may be large, even 
far from the spectrum [18] , and a better understanding of the resolvent properties of 
non-selfadjoint operators remains a challenging problem [29]. In this paper we are 
particularly interested in (polynomial) upper bounds in h for the resolvent. These 
bounds are especially useful when studying the stabilization problem, which deals 
with the rate of the energy decay of the solution of the damped wave equation on 
X : 

f (df + A g + a(x)d t )u(x,t) =0, (x,f)eXx(0,oo) 
[ ' \ u(x,0) = u a e H l (X), d t u(x,0) = ui£H°(X). 

It has been shown (see [53]) that if a > somewhere, then the energy of the waves, 

E(u,t) = - [ \d t u\ 2 + \Vu\ 2 dx 

satisfies E(u,t) *~ i ' 00 > for any initial data (uo,Ui) £ H 1 x H°. If some monotone 
decreasing function f(t) can be found such that 

E(u,t)^f(t)E(u,Q), 
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so-called strong stabilization occurs. It is not hard to show that this is equivalent 
to a uniform exponential decay : 3C, (3 > such that for any u solution of (|1.2j) . 

E{u,t) s$ Ce- 0t E(u,O). 

In pioneering works of Rauch, Taylor, Bardos and Lebeau it has been 

shown in various settings that strong stabilization is equivalent to the geometric 
control condition (GCC) : there exists T > such that from every point in E = 
= 1} c T*X, the bicharacteristic of P(h) reaches {a > 0} in time < T . By 
contrast, when the manifold X is no longer controlled by a, decay rate estimates 
usually involve additional regularity of the initial data. They take the form 

E(u,t)< f.(t)\Mw 

for s > and 

\\u\\^ = \\n(0)\\W + \\dtu(0)\\h- 
The question of exponential energy decay reduces to the study of high-frequency 
phenomena, in particular the issue of the spectral properties in the semiclassical 
limit h — > of certain non-selfadjoint operators approximately of the fornQ (jl.ll) . 
on a fixed energy layer. For instance, when geometric control holds, there exist 
h > 0, C, c > such that for z € [1 — 5, 1 + S] + i[— ch, ch] we have 

(1.3) \\Rz(h)\\ C (L^Li) <C/h, h<h Q . 

Standard arguments then show that this resolvent estimate implies the uniform ex- 
ponential decay for the energy. Similar arguments will apply in the case considered 
here, of resolvent estimates with loss. 

1.1. Motivation. While our main motivation for studying resolvent estimates for 
P{h) come from the stabilization problem, our approach in this paper is oriented by 
geometric considerations, as we explain now. As discussed above, in the presence 
of geometric control, the resolvent (P(h) — z)^ 1 enjoys a polynomial bound in a 
neighbourhood of size ch around the real axis, and this property implies exponential 
damping. When geometric control no longer holds, it is then a natural question 
to ask what type of estimate can be satisfied by \\(P(h) — and, crucially, in 

what type of complex neighbourhood of the real axis a resolvent estimate can be 
obtained. The properties of the undamped set 

(1.4) Af = {p E S*X : Vt 6 R, o o c tHp (p) = 0} 

are of central importance for this question. Here, H p denotes the Hamiltonian 
vector field generated by the principal symbol p = (jh{P(h)) = |£|g of the operator 
P(h), and S*X = p _1 (l) denotes the unit cosphere bundle. We remind the reader 
that the flow generated by H p in S*X is simply the geodesic flow. 

We now review some known results in the case Af ^ 0. In [7], the case when 
Af is a single hyperbolic orbit is analyzed, and a polynomial resolvent estimate 
for || (P — ) 1 1| is shown in a h/\ log/i|-size neighbourhood of the real axis. As a 
consequence, the energy decay is sub-exponential : with the above notation, one 
can get f s (t) = e _,3sv/ *. It is known from recent work [3] that this decay is sharp. If 

'Strictly speaking, in order to apply resolvent bounds to the damped wave equation, we also 
need the imaginary part of the Schrodinger operator to be mildly z dependent, with 

R z {h) = (h 2 A g + ihay/z- z)' 1 ; 

this will be handled by perturbation (see Corollary |4.3l Section [6] as well as references [9l l21ll23l ). 
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the curvature of X is assumed to be negative, and if the relative size of the damping 
function a is sufficiently large, then the resolvent obeys a polynomial bound in a 
size h neighbourhood of the real axis, and as a result, exponential decay for regular 
initial data occurs [27][28]. Indeed, the hypotheses in [28] is much more general, 
requiring only undamped sets of small pressure; the hyperbolic geodesic is a special 
case. We note also that for constant negative curvature, the need for an arbitrarily 
large damping function a has been recently removed by Nonnenmachcr, by using 
different methods [2"4] . 

A natural question raised by the above remarks, is the following: to what extent 
does the geometry of the trapped set alone determine a type of decay? In other 
words, given a trapped geometry, what type of resolvent estimate do we expect, 
and in what complex neighbourhood of the real axis? This amounts in many cases 
to a potentially rather crude decay rate for the energy, as it only depends on the 
structure of AT and not on the global dynamics of geodesies passing through the 
damping; in certain cases, however, our results can be seen to be sharpQ 

Motivated by the "black box" approach of Burq-Zworski [4] (cf. earlier work 
of Sjbstrand-Zworski [31]) as well as recent work on the gluing of resolvent esti- 
mates by Datchcv-Vasy |15j , we give a recipe for taking information from resolvent 
estimates obtained for a noncompact problem in which the set Af consists of all 
trapped geodesies — those not escaping to infinity — and investigate what these esti- 
mates imply for the compact problem with damping. In practice, as recent results 
of Datchev-Vasy [15] have shown the resolvent estimates on manifolds with, say, 
asymptotically Euclidean ends to be equivalent to estimates on a compact manifold 
with a complex absorbing potential substituting for the noncompact ends, we choose 
for the sake of brevity and elegance to take this latter model as our "noncompact" 
setting. As will be discussed below, these complex absorbing potentials have the 
effect of annihilating scmiclassical wavefront set along geodesies passing through 
them; this is why they are roughly interchangeable with noncompact ends, along 
which energy can flow off to infinity never to return. 

We thus formulate our main question as follows. Assuming that (X, g) and a are 
given, we consider a model operator of the form 

Px(h) = h 2 A g +iWi 

in which the damping is replaced by a complex absorbingpotential. We assume that 
this model operator enjoys a given "resolvent" estimate^ on the real axis: 

(1.5) || (PxC/i)-^- 1 1| ^C^P-, ze[l-S,l + S\. 

h 

Various examples of such estimates already appear in the literature, see for instance 
[5]IT]I5 ]I11|[T2 ,25 34,31)] and the references therein. Given (|1.5[) we then aim to obtain 
analogous estimates for the inverse of the operator with damping, i.e., on 

{P{h) - z)- 1 = (h 2 A g +iha- z)' 1 , 

when the complex absorbing potential has crucially been replaced by an 0(h) damp- 
ing term. In this paper, we address this question using a control theory argument 

t A very natural further question would then be : given a trapped geometry, what kind of global 
assumptions on the manifold can improve the crude decay rate obtained when only M is known? 

* We refer to this as a resolvent estimate owing to its close relationship with the estimate for 
the resolvent in scattering problems. 
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motivated by [4], together with a recently improved estimate on resolvents trun- 
cated away from the trapped set on one side [17] . and show that we obtain the 
same order of estimate as for the model operator. In the next subsection we state 
the precise results. 

1.2. Results. As above, we take H p to be the Hamilton vector field of p = |£| — 1 
and e tHp its bicharacteristic flow inside p~ 1 (l) = S*X (i.e., geodesic flow). We 
continue to take 

(1.6) JV = {p G S*X : Vi e R, a o e tH " (p) = 0}, 

and will add the further assumption that 



7r(A0 n suppa = 0, 

where n is projection T* X — > X. Thus there exists a non-empty open set 0\ such 
that supp(a) C X\0±, and 7t(jV) <e 0\. The following Theorem is our main "black 
box" spectral estimate. 

Theorem 1.1. Assume that for some 5 € (0,1) fixed and K € Z, there is a function 
1 s$ a(/i) = 0{h~ K ) such that 

\\{h?& g + xa-z)- r \\ L ^ L 2 

/or z€ [1 — <5, 1 + 5] . T/ien f/iere exists C,cq > smc/i i/iai 

IK^A^ + i/ia-z)- 1 !!^, s^C^, 

for z e [1-5,1 + 5]+ i[-co,c }h/a(h). 

When A/" = 0, one has a(/i) = 0(1), while for A/" ^ 0, one has a(/i) — >• oo as 
/i — > 0. As a general heuristic, the "larger" the trapped set is, the larger is a(h) 
when h — > 0, and the weaker the above global estimate is — see section [5] below for 
examples. 

Remark 1.2. As discussed above, instead of the assumption on the model operator 
h 2 A g + iWi with complex absorption, wc could just as well, by the results of |15j . 
have made an assumption on a model operator in which the set 0\ is "glued" to 
non-compact ends of various forms. In particular, it would suffice to know the 
cut-off resolvent estimate on the real axis for the limiting resolvent 

\\ x (h A -z + iQ) x\\ c{L 2 X 2 } < — 

for a localizer \ equal to 1 on 0\ and for A' the Laplacian on a manifold with 
Euclidean ends whose trapped set is contained in a set 0[ isometric to 0%. 

Remark 1.3. The hypotheses of Theorem 1 1 . 1 1 can be weakened to phase space hy- 
potheses, with a a pseudodifferential operator (as in [3D]). We have chosen to keep 
the damping as a function on the base in accordance with tradition and for the sake 
of brevity. 

Remark 1.4. The assumption that a(h) = 0(h~ K ) is of a technical nature. It does 
not appear to be too restrictive, however, since every known estimate for weakly 
unstable trapping satisfies this assumption (see Section [5] below). Indeed, if the 
undamped set N is at least weakly semi-unstable, the results of [TTJ suggest that 
in fact a{h) is always (^(/i^ 1 ^ 6 ) for any e > 0. 
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Remark 1.5. If a(h) is not of polynomial nature, the proof of Theorem 11.11 has to 
be slightly modified (see below). As a result, the final estimate we can obtain is 
weaker : there exists C, Co > such that 

\\(h 2 A a + iha - z)- 1 ^^ ^ C?-jp-, 

for z £ [1 - S : 1 + S] + i[-co, c ]h/a 2 {h). 

In section [5] wc describe three different settings in which our gluing results apply, 
in which the dynamics in a neighborhood of the trapped set are respectively 

(1) Normally hyperbolic 

(2) Degenerate hyperbolic 

(3) Hyperbolic with a condition on topological pressure. 

In addition to proving resolvent estimates in these settings, we discuss applications 
to decay rates for solutions to the damped wave equation. 

2. Operators with complex absorbing potentials 

In this section, we collect some standard results about operators with complex 
absorbing potentials. Such a potential has a much stronger effect than damping, 
namely (in the microlocal absence of forcing), that of annihilating semiclassical 
wavefront set completely along bicharacteristics passing through it in the forward 
direction. 

We collect basic results about the resolvent of an operator with complex absorb- 
ing potential. This includes both existence of the family and the basic propagation 
estimates, which tell us that the complex absorbing potential kills off wavefront set 
under forward propagation. We begin with the definition of the "resolvent:" 

Lemma 2.1. Let W £ C°°(X), W ^ and W not identically zero. Suppose that 

Pi = h 2 Ag +iW 

on X. Then (Pi — z)~ l is a meromorphic family of bounded operators on L 2 for 
all z £ C, analytic in the closed lower half-plane. 

Proof. We simply remark that h 2 A+l is invertible, and apply the analytic Fredholm 
theorem to conclude that there exists a meromorphic resolvent family. By the 
Fredholm alternative, any pole has to correspond to nullspace of Pi — z. Since 

Im (Piu, u) = (Wu, u) — (Imz)||u|| 2 , 

for u to be in the nullspace would imply that Imz 0; equality would further 
require that (Wu, u) = which is forbidden by unique continuation. □ 

Finally, we recall the microlocal bound of propagation through trapping by 
Datchev-Vasy [16] . as well as basic backward propagation of singularities in the 
presence of complex absorption (see also Lemma A. 2 of [25] for the latter). In this 
context we will say that a bicharacteristic 7 (by bicharacteristic we always mean an 
integral curve of H^ cpi in the characteristic set of Repi — Re z where p\ is the semi- 
classical principal symbol of Pi), or a point on 7, is non-trapped if W(j(T)) > 
for some T £ R, and is trapped otherwise. We say it is forward non-trapped if 
W("f(T)) > for some T > 0. In the terminology of |16j we say thaf| the resolvent 



§ We have opposite signs for imaginary parts of pi relative to 1161 . so incoming and outgoing 
are interchanged. 
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(Pi — z) 1 is semiclassically incoming with a loss of h 1 provided that whenever 
q £ X*X is on a forward non-trapped bicharacteristic 7 of Pi and / = 0(1) on 
7 | [0jT] then (Px - z)" 1 / is Oih- 1 ) at g. 

Lemma 2.2. (See EE!-) Suppose that P x = h 2 A g + i W on X, W > 0, and z G C 
smc/i i/iai Imz = 0(h°°), Re z E [1 — 5,1 + 5] for < <5 < 1 /iied. Assume that the 
resolvent is polynomially bounded, i.e., 

|| (Pi — z )~ 1 ||£( i 2 £2) ^ Ch~ K for some K. 

Then (Pi — z) _1 is semiclassically incoming with a loss of h . In particular, if 
W("f(T)) > for some T > and WFr(/) is disjoint from 7|[o,t]> ^ien (Pi — 
z)- 1 / = 0(/i°°) af g. 

Further, if \ € C°°(^) and T* uppx X contains no trapped points, then x{Pi ~ 

Z)-! X ZSO^" 1 ). 

A further result that is of crucial importance in avoiding losses in our estimates 
is the following result of [17] : 

Lemma 2.3. (See \PTj.) With the notation of Lemma \2.H\ if 

IkPi-z)- 1 !! <^ 

II ^ Z > \\C(L 2 .L 2 ) /j 

and if x € C°°(X) and T* contains no trapped points, then for some C > 0, 

ll(*- *>-'*! W^c^S. 

3. Propagation and damping estimates 
We now switch from complex absorbing potentials back to damping: set 



P = P(/i) = /i 2 A g +iha 



and consider the equation 

(P-z)u = f, z€ [1-5,1 + 5], 

We also set £ = p _1 (0) C T*X where p = 07,. (P). Let us start by recalling 
a classical result about propagation estimates (see, e.g. Theorem 12.5 of [37] for 
a proof by conjugation to normal form; an alternative is the usual commutator 
argument as described in [22] in the homogeneous setting and [8] in the semiclassical 
setting) : 

Lemma 3.1. Suppose q € £ and for some T > £/ie forward bicharacteristic 
exp([0, T]H p )(q) is disjoint from a compact set K. Then there are Q,Q' which 
are elliptic at q, resp. exp(TH p )(g), with WF^(Q') l~l K = such that \\Qu\\ ^ 

||g' M |j+c^ 1 ||/||,u = (p-z)- 1 /. 

4. Proof of Theorem 11.11 
For the moment, we consider only the case where z G [1 — 5, 1 + 5]. We have 

P = h 2 A g +iha 



and we additionally write 

Pi = h 2 A g 
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for the operator with damping replaced by absorption. Choose an open set Vi such 
that Af <s Vi <s 0\. Let Bi,(p <E C§°(X) be smooth functions with Pi|vi = 1, 
suppPi C 0\, (p = 1 on supp VPi and suppt^ PI Af = 0. We observe that Af 
satisfies the assumptions of Lemma 12.31 so that 

|| (Pi - «) -1 H| < Ch- x ^/^{K)\\ipu\\. 
Then, noticing that a and Pi have disjoint supports, we have 
\\B lU \\ = \\(P 1 -z)- 1 (P 1 -z)B 1 u\\ 
= || (Pi -z)-\P -z)B 1 u\\ 

( 4.i) = ii(p 1 -^)- i (p.i(p-^) + [p,pi]>ii 

(KPj - z)- 1 Pi(P - z)u\\ + || (Pi - z)~ V[P, Pi]HI 

< ^||s 1 (p-zH| + c A /^)||HI 

Since for p ^ A/" the curve e* Hp (p) passes through {a > 0}, each such bicharac- 
teristic curve must certainly enter the compact set X\0\. Thus by compactness, 
there exists eo > such that every such curve passes through {a ^ eo}. We now 
take a cutoff function x ^ with suppx C {a > £o/2}, and x = 1 whenever a ^ eo; 
hence every controlled geodesic passes through {x = !}■ 

We next recall a classical lemma concerning the propagation of singularities in 
the presence of geometric control. (See [3], which builds on a semiclassical version 
of [32], proved in [32].) This is a slight variation on Lemma \'3.1l and can of course 
also be proved using the original positive commutator argument (see [8]). 

Lemma 4.1. (See ^) Let U be an open neighbourhood of AT, x & C°°(X) as 
above. If Bq = (X) is such that WF^(Po) C T* X \ U , then for z real near 1, 

\\b u\\ < ^-\\(p-z)u\\ + \\xu\\ + o(h°°)\\u\\ 

Since supp<^ and supp(f — Pi) lie inside the controlled region, we can write : 
||(7 -PiHI s$ Ch^WiP - z)u\\ + C|M| +O(h 00 )\\u\\ 

and 

||HI Ch- l \\{P - z)u\\ + C\\ X u\\ + 0(h°°)\\u\\ . 

But 

= Ch~ 1 Im{(P-z)u,u) 
s£ Ch^WiP - z)u\\\\u\\. 

Starting with (|4.1[) . we deduce from the above inequalities that 
||Pi«|| 2 ^c(^-\\{P-z)uf + Ga{h)\\ipuf) 

< c {^rW( p - Z H 2 + ^(h)\\ X u\\ 2 + o(hn\\u\\ 2 ) . 
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Hence, we have 

\\u\\ 2 ^CQlBruf + Ul-B^uW 2 ) 



'o^(/i), 
h- 

{^\\(P-*)u\\ 2 + ^\\(P-z)u\\\\u\\ 



< c(^y\\(P-z)u\\ 2 + a(h)\\ X uf + 0(A°°)||«|| a 



■0(h°°)\\u 



<c(-y\\(P-z)u\\ 2 + — j^\\(P-z)u\\ 2 + e\\u\\ 2 

+ O(h 00 )\\u\\ 2> 



If e is small, we can absorb the last two terms in the above inequality on the left 
hand side, and we obtain 

M<C ^\\(P-z)u\\. 

Now simply observe that by the triangle inequality this bound is still valid if we 
add to z an imaginary part that satisfies 

|Imz| sC hoT l (h)C' 

for C such that C'C < 1, and this concludes the proof of the theorem. 

Remark 4.2. If a(h) is not of polynomial nature, then Lemma 12.31 cannot be used. 
As a result, the square root in Equation (|4.1[) must be removed. The rest of the 
argument is the same, and we end up with the estimate given in Remark 1 1.51 Note 
also that the energy decay rates for the damped wave equation are of course weaker 
than in the case where Lemma [231 can be applied. 



In order to apply Theorem 11.11 to the situation of the stationary damped wave 
operator, we now state a corollary where the Schrodinger operator depends mildly 
on z. 

Corollary 4.3. Let X be a compact manifold without boundary, and let P(h, z) be 
the modified operator 



P(h, z) = h 2 A g + i hy/za - 



z . 



Assume that for some 8 £ (0, 1) fixed there is a function 1 ^ a(h) = 0(h ), 
for some K € Z, such that 

\\ { h 2 A g + \a-z)-'\\ L ^ L ,^^, 

for z G [1 — 5, 1 + S] . Then there exists C,cq > such that 

\\P(h,z)-^ L2 ^C^, 

for 

z e [1 — Co/a(h), 1 + c /a(h)] + i[—c ,c ]h/a(h). 
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Proof. The real part of the perturbation is manifestly bounded by a small multiple 
of h/a(h) and can thus be perturbed away by Neumann series. It thus only remains 
to check that the size of the imaginary part of the perturbation: 

Rc(/ia — \f~zha) — Im z 

can also be made much smaller than h/a(h). 

Takc x /i= l+r+i/3 for r,(3 to be determined. Then z = (l+r) 2 -/3 2 +2(l+r) i/3. 
Then 

Re(ha - y/zha) - Imz = ha(l - (1 + r)) - 2(1 + r)/3 = eh/a(h) 
for e > small if, say, |/3| ^ eh/2a(h) and \r\ ^ Squaring, we obtain 

z e [1 — co/a(h), 1 + co/a(ft)] + i[— c , c ]/i/a(/i). 

□ 

5. Examples 

In this section, we briefly outline some known microlocal resolvent estimates, 
state the two different stationary damped wave operator estimates, and then draw 
conclusions about solutions to the damped wave equation (|1.2I) . 

5.1. A normally hyperbolic trapped set. In this section, we treat the case in 
which the trapped set is a smooth manifold in S*X around which the dynamics 
is normally hyperbolic. In this case, estimates for the resolvent with a complex 
absorbing potential have been obtained by Wunsch-Zworski [36] . A particular case 
of interest is the "photon sphere" for the Kerr black hole geometry, where the phase 
space is 6-dimensional A/" is a symplectic submanifold, diffcomorphic to T*S 2 — see 
section 2 of [35] for details on this application, and [33J for placing it in actual 
Kerr-de Sitter space. Another special case is of course that of a single hyperbolic 
closed geodesic (discussed further in the following section). 

The precise formulation of normal hypcrbolicity used here is as follows: we define 
the backward-forward trapped sets by 

r± = {p e T*X : Vt ^ 0,aoe tHp (p) = 0} 

Then of course 

A" = r+nr_, 

where we have now ceased to restrict to a single energy surface (so M C T* X is 
a homogeneous subset in view of the homogeneity of p) in order to employ the 
terminology of symplectic geometry more easily. 

We make the following assumptions on this intersection: 

(1) r± are codimension-one smooth manifolds intersecting transversely at Af. 
(It is not difficult to verify that T± must then be coisotropic and A^ sym- 
plectic.) 

(2) The flow is hyperbolic in the normal directions to K within the energy 
surface S*X: there exist subbundlcs E ± of Tjv(r±) such that at p £ S*X 

7Vr± =TN®E ± 1 

where 

d(exp(H p ) : E ± -> E ± 
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{a > 0} 




Figure 1 . The "peanut of rotation" . 



and there exists 9 > such that 

(5.1) ||d(exp(Hp)(u)|| < Ce^^lMI for all v £ £ T , ±t > 0. 

As discussed in [35], these hypotheses as stated are not structurally stable, but they 
do follow (at least up to loss of derivatives) from the more stringent hypothesis that 
the dynamics be r -normally-hyperbolic for every r in the sense of |20i Definition 
4]. This implication, and the structural stability of the hypotheses, follows from a 
deep theorem of Hirsch-Pugh-Shub [20] and Fenichel [P9] . 

As a consequence of the estimates of [36] for resolvents, we then obtain the 
following estimate for the damped operator: 

Theorem 5.1. Let (X,g) satisfy the dynamical conditions enumerated above. Then 
we have 

1 1 (ft 2 A 9 + i ha - z) - 1 1 1 L ^ L * < , 

for z £ [1 — 6,1 + 5] + i[-co, c ]h/\logh\. 

This estimate, or more precisely its refinement in Corollary 14.31 provide a corre- 
sponding energy decay estimate for solutions to the damped wave equation (|1.2p . 
In order to avoid irritating issues of projecting away from constant subspaces, etc., 
we assume that u = u(x, 0) = 0. 

Corollary 5.2. Assume the hypotheses of Theorem \5.1\ hold. and letu be a solution 
to (|1.2p with uo = 0, and u\ £ H s for some s £ (0, 2]. Then there exists a constant 
C = C s > such that 

E(u,t) ^C7 e - stl/2 / c || Ul ||^. 

A simple situation in which the hypotheses of Theorem 11.11 are satisfied is that 
of a connected compact manifold of the form X = Xq U X\ with X\ open and Xq 
isometric to a warped product R u x S@ ~ l with metric 

g = du 2 + cosh 2 u d8 2 . 

We take a £ C°°(X) to be identically 1 on X\ as well as equal to 1 for |x| > 1 in Xq. 
This class of manifolds thus includes the "peanut of rotation" shown in Figure Q] as 
well as its higher dimensional generalizations. 

Then the trapped set is easily seen to be M — {u = 0, £ = 0} where £ is the 
cotangant variable dual to u, and the function x = 2 — u 2 satisfies the convexity 
hypotheses. The flow on M is normally hyperbolic, with the stable and unstable 
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manifolds being the two components of the set 




) = M 2 *,-, 



i.e., by the intersection of the condition that energy and angular momentum match 
their values on J\f '. This is an example of a normally hyperbolic trapped set, and 
hence both parts of Theorem 15.11 apply. 

5.2. A trapped set with degenerate hyperbolicity. In this section, we study 
a variant of the normally hyperbolic case, in which the intersection of stable and 
unstable manifolds is no longer transverse, hence the results of |36j no longer apply. 
This is the case of a surface of rotation with a degenerate hyperbolic closed orbit. 

Our example manifold is a topological torus X = [—1, l] x x Sj, equipped with 
the metric 



and m is an integer ^ 1. This manifold has a "fatter" part and a "thinner" part. 
At the thickest, there is an elliptic geodesic, and at the thinnest part, where x = 0, 
there is an unstable geodesic, which we denote by 7. If m ^ 2, the Gaussian 
curvature is chosen to vanish to a finite order at the unstable geodesic, hence the 
geodesic is degenerately hyperbolic. If the Gaussian curvature is strictly negative 
(m = 1) in a neighbourhood of the thinnest part, the geodesic is non-degenerate; 
the geometry of a single closed hyperbolic geodesic has been extensively studied 
in [3J21GIJE1QII1E2H1] an d others. As is seen in the previous section, in this non- 
degenerate hyperbolic case the energy decays sub-exponentially with derivative loss; 
in [3], it is shown that the sub-exponential decay rate is sharp. Based on the sharp 
polynomial loss in local smoothing and resolvent estimates in [12] , Theorem 11.11 
shows that for the degenerate hyperbolic periodic geodesic, we have the following 
estimates. 

Theorem 5.3. Let X be as above, and suppose a(x) controls X geometrically 
outside a sufficiently small neighbourhood U ~D 7, so that M — {7}. Then 



As in the previous subsection, we deduce from this resolvent estimate an energy 
decay estimate for solutions to the damped wave equation. 

Corollary 5.4. Assume the hypotheses of Theorem \5. 3\ hold, and letu be a solution 
to p. 21) with uq = 0, and u\ G H s for some s £ (0, 2]. Then there exists a constant 
C = C s > such that 



(5.2) 



ds 2 = dx 2 + A 2 {x)d9 2 



where near x = 0, 



A(x) = (1 + \x\ 2m )^ 



IK^Ag + i/ia-z)- 1 ^! ==; c?i- 2m/(m+1) |M| 

for z G [1 - S, 1 + 8] + i[-co, c ]/i 2m /(" i+1 ) . 
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5.3. Hyperbolic trapped set with small topological pressure. In this sec- 
tion, we assume that the trapped set Af has a hyperbolic structure, and that the 
topological pressure of half the unstable Jacobian on the trapped set is negative. 
Roughly, this means that the set Af is rather thin, or filamentary: in dimension 
2, this is for instance equivalent to require that Af has Hausdorff dimension < 2. 
The simplest case to have in mind is a single, closed hyperbolic orbit. We then 
build on |25| to get resolvent estimates near the trapped set, which we extend to 
the global manifold with our different methods. 

We briefly recall here the above dynamical notions. By definition, the hyperbol- 
icity of Af C S*X means that for any p € Af, the tangent space T p Af splits into 
flow, stable and unstable subspaces 

TpAf = RHp ®E s p ®E u p . 

If X is of dimension d, the spaces E p and E p arc d — 1 dimensional, and are 
preserved under the flow map: 

Vt £ R, de tH *(E s p ) = E s ctHp(py de m »(E u p ) = £> p(p) . 

Moreover, there exist C, A > such that 

i) ||de tH *(v)|| < Cc- Xt |H|, for all v £ E s p , t.^0 
(5.3) it) \\de- tH »{v)\\ ^ Cc- Xt \\v\\, for aU v € E u p , t> 0. 

One can show that there exist a metric on T* X call the adapted metric, for which 
one can take C = 1 in the preceding equations. 

The above properties allow us to define the unstable Jacobian. The adapted 
metric on T* X induces a volume form O p on any d dimensional subspace of T(T*X). 
Using dp, we can define the unstable Jacobian at p for time t. Let us define the 
weak-stable and weak-unstable subspaces at p by 

Ef = E s p ® RHp , E u / = E u p ® RHp. 

We set 

Jn.) = det d e-^| = n^,A- A -A^^) ;= 

where (vi,. . . ,Vd) can be any basis of E^ Hp . . While we do not necessarily have 

J u {p) < 1, it is true that J"(/o) decays exponentially as t — > +oo. 

We denote by Pr^ the topological pressure functional on the closed, invariant 
set Af. We briefly recall a definition, see [3S], [35] for more material. If / is a 
continuous function on Af, n an integer and e > 0, define 

!n-l 
^exp^/oe fcH -(p) 
P es k=o 

where the supremum is taken over all the (e, n) separated subsets S. The topological 
pressure of / on Af is then 

Pw(/) := lim limsup - log Z n>e (f) . 
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Our main assumption here is that the topological pressure of \ log J u on M is 
negative, namely: 

Prolog J u ) <0. 
For some 6 > small enough, this imply the following resolvent estimate: 

(5.4) Vze [1-6,1 + 6], \\(h 2 A g +ia- z)-^ ^ <7^jp 

This estimate is already contained in |25j . modulo two minor simplifications 
in our case: the manifold is compact, and infinity is replaced with the absorbing 
potential xa, which control everything outside the trapped set - |15j shows explicitly 
that the estimate, with the more complicated geometry at infinity, of [2 5) implies 
the slightly simpler complex absorption result. Using Theorem 1 we immediately 
deduce the following result: 

Theorem 5.5. Let X be a compact manifold, and suppose that a ^ controls X 
except on N , which is assumed to be hyperbolic with the property that 

PrA^logJ 11 ) <0. 

For 6 > small enough, there is ho and Cq > such that for h ^ ho and z G 
[1-6,1 + 6]+ i[-c , c ] we have 

,2 A ,.. ^-Hi ^ -N^l 



\(h 2 A g +iha- z)- L \\ C- 



h 

In particular, there is no spectrum near the real axis in a region of size h/\\ogh\. 
As the resolvent estimate is the same order as that in Theorem 15. 11 we deduce the 
same energy decay estimates as in Corollarv l5.2l 



6. From resolvent estimates to the damped wave equation and 

energy decay 

In this section, we show how to move from a high energy resolvent estimate to 
an energy decay estimate for the damped wave equation, proving Corollaries 15.21 
and 15.41 To estimate the energy decay for the damped wave equation, as usual we 
rewrite it as a first-order evolution problem : if it — (u, dtu) one can write (|1.2[) as 



-ild 
i a 



(6.1) d t u = iBu, B = 

The evolution group e ltB maps initial data (u ,ui) e H := i? 1 (X) x H°(X) to a 
solution (u,dtu) of ()6.1|) where u solves (|1.2|) . For s > 0, define 

IMU := \\uo\\m+°(X) + \\ui\\h°(x) 
It is not hard to see that if we can prove 

( 6 - 2 ) IK^-^rll^^po^/c*) 

then we can deduce a decay rate for the energy : 

E(u,t) ^ f(t)\\u\\ 2 s 

It turns out that we can obtain bounds such as (|6.2|) if we have estimates on the 
high-frequency resolvent (A — |A| — > oo. 
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0=<AeC:|imA|< 



To see this, we recall the following setup from [9]. Now suppose (A — 23) 1 
continues holomorphically to a neighbourhood of the region 

Ci, | Re A| < C 2 

P(|ReA|), |ReA|^C 2 , 

where P(|ReA|) > and is monotone decreasing (or constant) as |ReA| — > oo, 
P{C2) = Cij and assume for simplicity that 3f2 is smooth. Assume 

(6.3) ||(A-23)- 1 || if ^ <G(|ReA|) 

for A G O, where G(| Re A|) = 0{\ Re \\ N ) for some N > 0. 

Theorem 6.1. Suppose 23 satisfies all the assumptions above, and let k € hi, fc > 
A/ + 1. TTien /or any f (t) > 0, monotone increasing, satisfying 

(6.4) F(i)( fc+1 )/ 2 sC cxp(tP(F(t))), 
i/iere is a constant C > suc/i £/ia£ 

(6.5) 



«S CF(t)- k/2 . 



(l-i23) fc 

In all cases considered in this paper, we have semiclassical resolvent estimates 
\\(h 2 A g + i^ha- z)-^^^^ < z ~l + ih/a(h), 

If we rescale 

then our resolvent estimates become 

ii , * . o> 1 ,, adr! -1 ) 

IKAa+iro-T 2 )- 1 !!^^ < 7 - 

for Imr ~ (a(| Rcr| -1 )) -1 . By interpolation, this implies for ^ j ^ 2, 

IKA^ + ira-T 2 )" 1 !!^^ < Irr^drl" 1 ). 

Hence, with 23 as above and H = H 1 x a simple calculation yields 

||(A-B)- 1 || H _ >H <a(|A|- 1 ). 

For CorollaryEH we take aQAl" 1 ) = log(2 + |A|). Then jfe = 2 suffices, P(r) = 
log - (r), and hence we may take 

F(t)=c tl/2 / c . 

This recovers the endpoint estimate s = 2. To get the intermediate estimates for 
s £ (0, 2) we interpolate with the trivial estimate 

E(u,t) < P(w,0). 

For CorollaryEH we have a(|A| _1 ) = ^i™- 1 )/^ 1 ) , N = (m-l)/(m + l) < 1, 
so that k = 2, and P(r) = r (i-™)/(™+i). We try 

F(t) " E?W' 

and insist 

t 3.s/2 log -3 9 /2 (i) ^ exp (^(l— )/(m+l) logg ( m -l)/(m+l) W)) 
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which is satisfied if 

m + 1 
S= (m=i) 

and 

_ 3(m + l) 2 
9 ~ 2(m- l) 2 ' 

Again interpolating with the trivial estimate proves the Corollary. 
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